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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to the processing of positioning signals, such as 
10 signals received from the NAVSTAR Global Positioning System (GPS). In particular, 
the present invention relates to detecting positioning signals, by first using a less 
sensitive, more robust technique, followed with further processing using a more sensitive 
technique. 

2. Discussion of the Related Art 



15 Algorithms for the detecting positioning signals, such as GPS signals, differ in 

their sensitivity, robustness, ease of implementation, and computational resources 
requirements. Often, sensitive algorithms are computationally more demanding and less 
robust than less sensitive algorithms. 

Two well-known algorithms for signal detection are coherent integration and non- 
20 coherent integration. In coherent integration, the received positioning signal is 

multiplied variously with locally generated replicas of the expected satellite signal of 
different code phases and carrier frequencies, and each product is integrated for a time 
period equal to the signal capture time. Following this integration step, the magnitude of 
each integrated product is computed and, after further processing that may include 
25 interpolation, the magnitudes are each compared to a threshold to ascertain whether 
the positioning signal to be detected is present. Non-coherent integration is similar 
to coherent integration, except that the magnitudes or square magnitudes are 
computed at regular subintervals of the signal capture and the results are then 
summed, processed, and finally compared to a threshold value. 

30 In both the coherent and the non-coherent integration algorithms, multiple 

trial code phases and frequencies are tested. In coherent integration, the trial 
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frequencies (the "frequency bins") are spaced by — , where Tis the interval of 

coherent integration. Compared to non-coherent integration, coherent integration 
has a higher sensitivity, at the expense of a lesser robustness and a higher 
computational requirement. Coherent integration methods have been developed that 
5 trade robustness for less computational complexity. Coherent integration may also 
search over additional parameters such as oscillator phase acceleration. 

In certain environments, such as in urban locations, the positioning signals 
may be attenuated because of propagation through walls, reflections off of buildings 
and other signal distortions. Typically, because of the spatial distribution of the 

10 positioning signal sources, the positioning signals from each of the sources will be 
subject to different amounts of attenuation and reflection. As a result, it is rare that 
two positioning signals are attenuated by the same amount. Indeed, the signal 
powers of the received sources typically have a distribution in which there is 
considerable difference between the strongest positioning signal and the other 

15 signals. 

Numerous techniques for GPS signal detection are disclosed in the following 
U.S. patent applications: 



20 



1 . Signal Acquisition using Data Bit Information (Serial Number 
09/888,228 filed June 22, 2001, hereinafter referred to as "Application 
228"), 



2. Synthesizing Coherent Correlation Sums at One or Multiple 
Carrier Frequencies Using Correlation Sums Calculated at a Coarse Set 
of Frequencies (Serial Number 09/888,227 filed June 22, 2001, hereafter 
referred to as "Application 227"), 



25 



3. Extracting Fine-Tuned Estimates from Correlation Functions 
Evaluated at Limited Number of Values (Serial Number 09/888,338 filed 
June 22, 2001, hereafter referred to as "Application 338"), 



30 



4. Determining the Spatio-Temporal and Kinematic Parameters 
of a Signal Receiver and its Clock by Information Fusion (Serial Number 
09/888,229 filed June 22, 2001, hereafter referred to as "Application 229"), 



5. 



Determining Location Information Using Sampled Data 
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Containing LocationDetermining Signals And Noise (Serial Number 
09/888,337 filed June 22, 2001, hereafter referred to as "Application 337"), 

6. Method for optimal search scheduling in satellite 
acquisition (Serial Number 10/126,853, filed on April 19, 2002, hereafter 
referred to as "Application 853"), 

7. System and method to estimate the location of a receiver 
in a multi path environment (Serial Number 10/237,556, filed on 
September 6, 2002, hereafter referred to as "Application 556"), 

8. System and method to estimate the location of a receiver 
(Serial Number 10/237,557 filed on September 6, 2002, hereafter referred to 
as "Application 557"), and 

9. Multifunction device with positioning system and shared 
processor (Serial Number 10/286,360 filed on November 1, 2002, hereafter 
referred to as "Application 360"). 

The above patent applications are incorporated by reference herein in their 
entireties. 

SUMMARY OF THE INVENTION 

The present invention provides a systematic method for acquiring positioning 
signals, such as global positioning system (GPS) signals, using different signal 
detection algorithms at different stages of signal detection. According to one 
embodiment of the present invention, a method for detecting multiple positioning 
signals includes first detecting a first positioning signal using a robust but less 
sensitive signal detection method, such as non-coherent integration. Based on the 
signal parameter values that allow detection of the first positional signal, detecting a 
second positioning signal using a more sensitive method, such as coherent 
integration. In this manner, by capturing the strongest signal first using a robust 
method, signal detection parameter values can be narrowed to allow a subsequent 
signal acquisition using a more sensitive, albeit more computationally intensive, 
method. 

In one embodiment, the values of one or more signal detection parameters can 
be determined upon acquisition of the first signal. Such signal detection parameter 
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determination includes aligning the first positioning signal relative to a navigation 
message embedded in the first positioning signal, determining an oscillator 
frequency, and determining an oscillator phase acceleration. 

Non-coherent integration for capturing the initial positioning signal can be 
5 achieved using correlation values computed based on predetermined-size blocks of a 
received signal, such as 4 millisecond blocks. A search for oscillator frequency or 
oscillator phase acceleration can be performed in the frequency domain. Prior to 
signal detection using a more sensitive method (e.g., coherent integration), the effects 
of any navigation message in the positioning signal can be removed. In one 
10 embodiment, the effects of a navigation message is removed using a method 
involving doubling an estimated phase angle. 

Unlike methods and systems for processing positioning signals which utilize a 
single detection algorithm and which yields a particular sensitivity and robustness for 
a given computational cost, the methods of the present invention allows combining 

15 multiple detection algorithms to exploit the variation in the signal powers of the 
received positioning signals that commonly exist in the urban environment. By 
searching initially using a less sensitive and more robust algorithm, and upon the 
detection of the first positioning signal, computing additional information that is 
common to all the received positioning signals using the first detection, the search for 

20 the remaining positioning signals can use a more sensitive algorithm with greatly 
reduced computational requirements. For a mobile GPS device, the savings in 
computation intensity can often be large enough to capture positioning signals in an 
environment that makes such capture unattainable for the capabilities of such a 
device. 

25 When the received signal powers have a sufficiently wide distribution, the 

present invention combines the best attributes of the various detection algorithms; 
namely, it enjoys almost the same sensitivity of the most sensitive algorithms and 
the robustness and computational efficiency of the less sensitive algorithms. In 
urban environments, signal powers exhibit a sufficiently wide distribution to make 

30 this invention highly desirable. 

The present invention is better understood upon consideration of the detailed 
description below and the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates, generally, a method according to one embodiment of the 
present invention. 

Figure 2 illustrates in further detail a method in according to one 
5 embodiment of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 illustrates, generally, a method in accordance with one embodiment 
of the present invention. As shown in Figure 1, an initial search at step 101 uses a 
relatively less sensitive algorithm to find the strongest positioning signal. At step 

10 102, the initial search is checked to determine whether or not signal detection is 
successful. If signal detection is not successful, signal detection terminates at step 
103. If signal detection is successful, at step 104, the strongest signal thus detected 
is used to compute one or more parameters common to all positioning signals (e.g., 
the oscillator frequency offset and phase acceleration 1 ). At step 105, the computed 

15 values of the common parameters are used in the search for the remaining signals 
using a more sensitive algorithm. 

Figure 2 illustrates in further detail a method according to one embodiment 
of the present invention. As shown in Figure 2, at step 201, an initial search for a 
positioning signal is performed using non-coherent integration. In this embodiment, 

20 a received GPS positioning signal is divided into 4 milliseconds (ms) blocks. In 
general, block sizes between 1 ms and 14 ms may be chosen; choosing a shorter 
block allows a less fine frequency search, while choosing a longer block yields 
better sensitivity. Such non-coherent integration does not require a precise 
navigation message alignment, and is robust to oscillator phase accelerations, and 

25 requires a spacing between trial frequencies in the order of 125 Hz. Typically, a 

non-coherent integration over one second of a GPS signal is about 6 dB less sensitive 
than a coherent integration over the same one second of the GPS signal. As the 
number of trial frequencies to search is proportional to the reciprocal of the 
integration period times the uncertainty in the code phase, if one uses 8 ms blocks 

30 for non-coherent processing, an improved sensitivity can be achieved, but such non- 
coherent processing would require either a better frequency estimate or testing 
additional trial frequencies. 

1 Phase acceleration is the rate of change of oscillator frequency 
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At step 202, the initial search is checked to determine if a positioning signal 
is detected. If a positioning signal is not detected, the search terminates at step 203. 
Alternatively, at step 204, upon successfully detecting an initial positioning signal, 
the detected positioning signal is used to estimate one or more signal parameters such 
5 as oscillator frequency, oscillator phase acceleration, and the navigation message 
alignment. If an approximate location of the receiver is known, an absolute time 
stamp can be calculated which can then be used to reduce the code phase 
uncertainties of the other positioning signals to be detected. 

According to one embodiment, at step 204, complex 1-ms segments of the 
10 detected signal are correlated with a reference signal of various trial carrier 

frequencies and code phases, as described in the copending patent applications 
incorporated by reference above. The 1-ms correlation length is selected as a 
convenient length for aligning the navigation message. Navigation message 
alignment can be estimated using, for example, the technique disclosed in U.S. Patent 
15 6,295,023, entitled, "Methods, mobile stations and systems for acquiring global 
positioning system timing information" by Scott Bloebaum, which is hereby 
incorporated by reference in its entirety. To align the navigation message, a first 
step determines a 20-ms navigation bit boundary. The navigation bit boundary 
positions can be determined using non-coherent integration using coherent blocks of 
20 20ms, and choosing the boundary position that maximizes the integration. 

Individual navigation bit positions can then be determined by testing different 
alignments of the navigation message, using coherent or non-coherent integration 
(with block sizes greater than 20 ms), and choosing the alignment that yields the 
highest integration value. 

25 In the presence of oscillator phase acceleration, a long coherent integration is 

difficult. Instead, for navigation message alignment, non-coherent integration using 
block sizes between 20-100 ms can be used or, alternatively, a coherent integration 
that also searches for phase acceleration can be used for navigation message 
alignment. If non-coherent integration is used, a block length which is a multiple of 

30 20 ms is preferably not used so as to avoid ambiguities. 

If navigation message alignment is not required or if a small sensitivity loss is 
acceptable, correlation lengths greater than 1 ms may be used to reduce memory 
requirements. The resulting complex correlations can then be used to refine the 
estimate of oscillator frequency, oscillator phase acceleration, the navigation 
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message alignment and, if an approximate location of the receiver is known, an 
approximate time-stamp. 

After alignment of the navigation message, oscillator phase acceleration can 
be estimated and estimates of oscillator frequency can be improved, using various 
trial oscillator phase accelerations and oscillator frequencies, and finding the values 
that provide the best coherent integration over the 1-ms correlations. In other words, 
if Sk denotes the 1-ms complex correlations for k =1, 2,3,.. .N, then the oscillator 
frequency / and phase acceleration a are given by the trial frequnency and the trial 
phase acceleration that maximize the following expression: 



10 max / a £ 



N 

2v(Jk+0.Sak 2 ) 



Note that this search for / and a operates on the 1-ms correlation values, and 

hence, for one second of GPS signal, each coherent integration involves 1000 

correlation values. Thus, even simple brute force methods are not computationally 

onerous. Even then, after aligning the navigation message, the 1-ms complex 

15 correlations can be combined to further reduce the number of values for integration. 
One such combination combines M 1-ms correlations into a single value t k using: 



2jg(((k-l)M+i)f +0.5a ((k-\)M+i) 2 ) 



t = Y <> e 2 



/-I 



where / and a are the current estimates of oscillator frequency and phase 
acceleration, respectively. 

20 Alternatively, instead of using a time domain method, the search for 

oscillator frequency and phase acceleration at step 204 may be carried out in the 
frequency domain. A Fast Fourier Transform (FFT) may be used to test many trial 
frequencies simultaneously. Trial phase acceleration can be tested in either the 
time domain by first multiplying each correlation by its phase acceleration or in the 

25 frequency domain by convolving the FFT of each correlation by the transform of 
the trial phase acceleration value. 

Navigation message alignment can also be performed after oscillator 
frequency and phase acceleration estimates. To estimate oscillator frequency and 
phase acceleration before aligning the navigation message, the effects of the 
30 navigation message is first removed. It is tempting to remove such effects by taking 
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the magnitude of the received signal. However, since taking magnitude of a signal 
destroys phase information, taking magnitude of a signal is therefore unsuitable for 
estimating phase acceleration. Since the correlation values are complex, an 
alternative method to eliminate the effects of the navigation message is to double the 
5 phase angle. For a given phase angle, doubling the phase angle is equivalent to first 
adding 180 degrees to the phase angle and then doubling the sum. Since adding 180 
degree is equivalent to reverse the sign of a complex number, doubling the phase 
angle removes the navigation message, without removing the phase angle 
information. Such a method can thus be used to estimate phase acceleration. 

At step 205, coherent integration can be used to search for additional and 
weaker positioning signals using, for example, methods disclosed in Applications 
557 and 227 and a longer coherent integration. The estimates of parameter values 
computed in step 204 are essential for the robustness and low computational cost of 
step 205. If an accurate timestamp is obtained at step 204, search for additional 
positioning signals can be carried out user a finer estimate of code phases. In one 
embodiment, coherent integrations of 2-10 ms (the integration interval can be 
determined based on the residual frequency uncertainty) are stored in a two- 
dimensional memory array in which the rows correspond code phases and the 
columns correspond to time samples. After correcting for the navigation message, a 
FFT of the data of each row allows the search to be carried out in the frequency 
domain. Further, refinements of code phase and frequency can also be performed 
using interpolation, as described in Application 338. 

By using non-coherent integration to detect the positional signal of the 
strongest satellite, the above algorithm has greatly reduced computational cost. 
25 From this strongest satellite, estimates for absolute time, oscillator frequency, and 
oscillator phase acceleration can be obtained. The time and oscillator frequency 
allows the search spaces of carrier frequency and code phases to be narrowed 
substantially. The estimates of the frequency and code phases thus obtained allow 
long coherent integration to facilitate searching positional signals from additional 
30 satellites. 

The above detailed description is provided to illustrate the specific 
embodiments of the present invention and is not intended to be limiting. Numerous 
variations and modifications within the scope of the present invention are possible. 
The present invention is set forth in the following claims. 
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